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Abstract 



' We determine the correlation between Polyakov loops in three dimensional 

■ SU(3) gauge theory in the confined region at finite temperature. For this 

purpose we perform lattice calculations for the number of steps in the 
K/i , temperature direction equal to six. This is expected to be in the scaling 

^ ■ region of the lattice theory. We compare the results to the bosonic string 



model. The agreement is very good for temperatures T < 0.7 Tc, where 
Tc is the critical temperature. In the region 0.7 Tc < T < Tc we enter the 
critical region, where the critical properties of the correlations are fixed by 
universality to be those of the two dimensional three state Potts model. 
Nevertheless, by calculating the critical lattice coupling, we show that the 
ratio of the critical temperature to the square root of the zero temperature 
string tension, where the latter is taken from the literature, remains very 
near to the string model prediction. 
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1 Introduction 



Three dimensional SU(3) gauge theory has many properties in common with the 
corresponding four dimensional theory. At zero temperature lattice calculations 
show the confinement of heavy quarks in a linear potential. Furthermore, they 
predict a mass gap and a nontrivial glueball spectrum. At finite temperature 
there is a deconfining phase transition. In contrast to SU(3) in four dimensions the 
transition is second order. Lattice calculations of the critical indices are consistent 
with the transition being in the universality class of the two dimensional three 
state Potts model p[], as expected from general arguments[2]. 

In a previous paper, we have studied in detail the equation of state in the high 
temperature phase by performing extensive lattice calculations, as is possible in 
three dimensions, and thus obtaining very precise results [3j. In this note we 
describe a study of the theory in the confined phase below the transition. We 
work at A^^ = 6, which from the experience gathered in [3] we believe is in the 
scaling region. Thus our result should be relevant for the continuum limit. In 
particular we calculate the correlations between Polyakov loops. It has been 
proposed long ago that these can be described by an effective string model, the 
Nambu-Goto bosonic string [H [5]. In this model the temperature dependent 
string tension becomes zero at a critical temperature corresponding to the 
Hagedorn temperature of the string theory. The critical index for the approach 
of the string tension to zero is v = 1/2. This index is different from that of 
the two dimensional three state Potts model, which has u = 5/6. In fact the 
above predictions from the string model are independent of the gauge group, but 
only depend on the dimensionality of space. Early investigations have shown that 
nevertheless the critical temperature is not far from the prediction of the Nambu- 
Goto model [H El El- It has also been shown that for low enough temperature 
the correlations between Polyakov loops are very well described by the effective 
string model [8]. Finally, field theory investigations have proven that the first 
three terms of an expansion in of the string model are in fact universal results 
for a fiuctuating bosonic string P |TT] . But the expansion of course does not tell 
anything about the existence of a singularity at T^. 

Here, from the lattice, we extend the earlier results on the determination of the 
ratio between the critical temperature and the square root of the string tension. 
Furthermore we study in detail the behaviour of the Polyakov loop correlations 
in the region Tc/2 < T < Tc. For this purpose, we have taken data for a dense 
set of temperatures, varying the temperature by changing the lattice coupling 
constant. This method is complementary to the one used in [8J, where the lattice 
extent in the temperature direction is varied. As will be shown, where our results 
overlap with those of 0, they agree with each other. It is also very interesting 
to study in detail the critical region, near Tc, where we do not expect the string 
model to describe the correlations. Instead, in fact, critical scaling is found to be 
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a very efficient tool to describe the correlation length between Polyakov loops. 
We postpone the discussion of these issues to a separate publication [T^ . 

In the next section we present the lattice setup and define the quantities of 
interest in the string context. The simulations are described in section 3, together 
with the determination of the critical temperature for a lattice extent in the 
temperature direction equal to 6. Section 4 is devoted to the comparison of the 
results obtained with the string picture, by reference to the Nambu-Goto model 
and to field theoretic approaches. A summary and conclusions are given in a last 
section. 



2 Polyakov loop correlations at finite tempera- 
ture. Lattice setup. 

We simulate the three dimensional SU (3) theory, regularized on a finite euclidean 
lattice with lattice spacing a, Nr points in the (inverse) temperature direction, 
defined as the direction, Ns points in the two space directions 1, 2 and periodic 
boundary conditions in all directions. We use the standard Wilson action: 

5(f/p)=/55^(l-^ReTrf/p), (1) 
p 

where P denotes one of the 3A^t- x Ng plaquettes on the lattice. Up is the product 

of the [/-matrices around the plaquette, and /? is the lattice coupling constant. 
From the classical limit of the lattice action we may write 

f^=—2 (2) 

where is the (dimensionful) gauge coupling constant of the continuum theory. 
We define the temperature and the volume of the lattice by 



^ = (3) 

V = {aNsf. (4) 

In the confined phase, we measure the correlations between Polyakov loops wind- 
ing around the temperature direction: 

L(xi, X2) = Tr JJ^ [/^ (f + ne^) , (5) 

n=0 
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where Ur denotes the hnk in the time direction, whose origin is located in space 
at X = (xi,X2) and eV is the unit vector in the time direction. In 2-dimensional 
space, we define the correlation function G{z) between two lines at distance z 
from each other by 



^^^) = ^[ J2 ^e{L{xi,X2)L*{xi + z,x'^)) + 

^ \xi,X2,x'2 

X (6) 

Re {L{xi, X2)L*{x[,X2 + z)) J 

Xl,x'^,X2 / 

Throughout the paper, the symbol (...) denotes an average over a set of gauge 
configurations. Due to periodic boundary conditions, any coordinate is under- 
stood modulo Ns, and all sums run over intervals of length Ns- Using the discrete 
lattice symmetries, we may rewrite G{z) as 

G{z) = ^Y.^e{mQ)L\z,X2)) (7) 



X2 



In the forthcoming analysis, it will be assumed that for z large enough (see details 
in section 3), contributions from excited states of the 2-dimensional system of 
Polyakov loops can be neglected, and G represented by the contribution of its 
ground state energy only. In this situation, we will use the parametrization 

G{z) = b cosh ( m(^ - z)) + c. (8) 
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Using ([3]), the ground state energy, here denoted m in lattice units, is related to 
its physical value M by 

m = Ma = ——. (9) 

In the right hand side of ([8]), 6 and c are assumed to be constants for given lattice 
parameters. The first term is the lattice expression of G for a free bosonic field, 
whose normalization determines b. In our region of investigation the constant c 
is essentially consistent with zero, and always negligible in practice. 

From now on, we will interpret M(T) as the ground state energy of a fiux 
tube, writing 

^f^ ^ ^ (10) 

where cr(T) is a temperature dependent string tension. We keep this definition 
of <j{T) for any T, and from Eqs. (Pl9l) rewrite the latter equation as 

^ = mN^- (11) 
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In order to study the temperature dependence of cr(T), and compare our data 
with the string picture it is convenient to use the zero temperature string tension 
(To as a scale, i.e. using Eqs. (piTT]) to write 

T 1 

(12) 



The quantity a^/oo has been recently measured from numerical simulations using 
large Ng lattices [HIB], leading to high accuracy results at /3 = 8.1489, 11.3711, 
14.7172, 18.131, 21 and 40. Because we need it for a dense set of (3 values in the 
10 to 22 range, we determined for any (3 by fits to the above numerical 

data. After some trials, we finally retained the following parametrization 

FM ^ -v^=|f^. (14) 
/i = 3.3257 z = 1.99 p = 3.69 (15) 

Not only all the data points are perfectly fitted, but it is interesting to notice that 
the zero and pole positions z and p in Fo-q suggest the existence of a cross-over from 
a weak to a strong coupling regime, for /5 of the order of a few units, well below 
the lowest (3 value needed. Moreover, we use this function only for interpolating 
accurate data. We do not quote the errors on the fitted parameters, as they are 
anyway highly correlated. What is important is the error in the interpolating 
function. We found that the absolute error on the inverse of the function Fo-g(/3) 
is approximately constant and equal to 0.002. 



3 Simulations and results 

The simulations were done in a standard way using the Wilson action ([T]). We 
used one heathbath sweep followed by four overrelaxation sweeps. The Polyakov 
loops correlations were measured according to formula ([6]). The measurements 
were performed every five combined (heathbath and overrelaxation) sweeps. At 
least 20,000 measurements were kept for every f3 and A^^ combination, their 
average number being of the order of 60,000. These measurements are not inde- 
pendent, and the autocorrelation time for the correlation functions was found to 
be 10 in the worst case {jS = 20), resulting in a minimum of 2,000 independent 
measurements (/5 = 20, Ng = 96). The ground state energies m were obtained by 
fitting formula ([H]), and their errors estimated using the blocked boostrap method 
as now described. 

The total set of measurements for given values of (3 and lattice size was first 
divided into blocks of 500 measurements, in order to take care of the autocor- 
relation time. Then 100 bootstrap samples were generated, each one resulting 
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from drawing a new set of blocks randomly from the original sample. For each 
bootstrap sample we calculated the average value of the correlation function. The 
average of those averages was taken as the final value of the correlation function, 
and their standard deviation as the corresponding error. 

The errors on the fitted parameters were obtained by repeating this procedure 
once more. For each bootstrap sample, we fitted formula ([8]) for z in an interval 
[zo,^s/2], using the value calculated from the errors found in the previous 
step. We so obtained 100 sets of fitted parameters m, 6, c and used the standard 
deviation between them as an estimate of their errors. 

Eq. ([8]) retains the contribution from the ground state only. In order to control 
the possible effect of higher states, we measured the dependence of the fitted 
masses on the zq parameter, and found that the systematic error associated with 
the choice of zq becomes smaller than the statistical error for zq > 4. In table [H 
we list the final values and statistical errors of the parameters. 

The value jSc of the lattice coupling at which the transition takes place has 
been estimated using the Binder cumulant approach [13], following the method 
described in [H] and applied to the variable |Lp, where 



72 ^(^l'^2)- (16) 



■3 Xl,X2 



Using this method, we obtain 

/3c = 21.36 ± O.Ol^'*-) ± (17) 

The systematic error quoted covers small differences observed if other methods 
are investigated and uncertainties due to unknown subleading terms in the scaling 
ansatz. These issues will be discussed in detail in the forthcoming paper on the 
critical region [12]. Here, the present systematic error does not affect our result 



on Tc/ y (t(0), which is the relevant quantity in the forthcoming comparison with 
the string model. Indeed, from the result quoted above and Eq. [121 we obtain 

T 

0.976(15), (18) 



V^(0) 

where most of the error comes from the error on the function F(j{(3). 



4 The string picture 

In a fundamental paper, Pisarski and Alvarez [3] long ago related the finite temper- 
ature behaviour of gauge theories in the confined region up to the phase transition 
to the finite temperature behaviour of a bosonic string randomly fluctuating in 
the transverse directions. If the string is supposed to be the Nambu-Goto string, 
the formula for the temperature dependent string tension is 
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= , ^ , (19) 

where D is the space-time dimension of the gauge theory. As we can see from the 
formula, the temperature dependent string tension is not dependent on the gauge 
group, but only on the number of transverse dimensions. The string tension goes 
to zero at a critical temperature 

T, = VSc^o/tt (20) 

which coincides with the Hagedorn temperature of the string model. The ap- 
proach to this singularity has the mean field behaviour with the exponent v = 1/2. 

Later, it was shown that in any dimension, the first term in the expansion in 
T'^/y/oQ is universal in any dimension[9], i.e. 

a(T) tt(D-2)T^ , , 

= 1- ^ ^ ^ +••• 21 
(Jo 6(Jo 

In three dimensions there is a special situation, as here the first three terms are 
universal and coincide with the developement of the expression for the Nambu- 
Goto string truncated to this order [TT] : 

= 1 h ... (22) 

(7o 6(Jo 72(Jo2 432(Jo3 ^ ' 

Of course, this polynomial expression contains no information about a singularity 
where (t(T) vanishes, such as that of the full Nambu-Goto expression. 

In figure [T] we show the quantity (t(T)/(To versus T / ^/oq as measured on the 
lattice using Eqs. fll2lll3|14l) . The data points correspond to various Ns values 
and cover the domain 0.5 < T / ^Joq < 1. We do not extend the domain to smaller 
T, because we may enter the strong coupling region. In the figure we have also 
included the two points from [S], which are in this region. As one can see our data 
agree with these points. From the figure one can also see that the NG-model is 
consistent with the lattice QCD results up to about T/ y/ao=0.7, above which the 
corresponding curve is definitely above the data. We also show the two curves 
corresponding to Eqs. fl21ll22p . both higher in the graph in the whole region (since 
any term of the full expansion of Eq. (1191) are negative). The second curve cannot 
be distinguished from the full NG model up to T/y/oo = 0.7. To fix the ideas, 
we just remark that if the last term in Eq. fl22|) were aboue 3 times larger, the 
agreement with the numerical results would be extended up to T/ctq about 0.8. 

We checked that the spatial lattice sizes Ns used in figured] were large enough 
for the above conclusions not being altered by any finite size effects for T/ ^/oq < 
0.9 . In practice, below T j ^/oq ~ 0.85, m is always larger than about 0.1 and the 
Ns values used at least 6 times the correlation length mT^ . In turn, at the largest 
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Figure 1: Finite temperature string tension in units of the zero temperature string 
tension as a function of temperature for various Ns values. Solid red squares 
come from [8]. The solid line represents the Nambu-Goto expression (1191) . and 
the dashed (dashed-dotted) line its expansion to the first (third) order in T^. 
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T 


Ns 


m 


b 


c 


12.5 


0.5251 


96 


0.5055(50) 


0.10(^2 ) 


e-ii 


-0.0000045(34) 


13 


0.5509 


64 


0.4540(27) 
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e- ( 


0.0000031(58) 


13.5 


0.5766 


64 


0.4029(30) 




e-6 


0.000005(11) 


14 


0.6023 


64 


0.3620(18) 


v.viyo ) 


e-6 


0.00001(1) 


14.5 


0.6279 


64 


0.3207(20) 


n OA Q(T A\ 

U.z4o(^i4 ) 


e-o 


0.000008(20) 


14.7172 


0.639 


64 


0.3068(14) 


n AHA ( T R\ 
U.4U4(^iO ) 


e-o 


0.000008(14) 


15 


0.6535 


64 


0.2866(12) 


U.oUo(^Zo ) 


e-o 


0.000000(16) 


15.5 


0.679 


64 


0.2564(14) 


0.234(10) 


e-4 


-0.000008(27) 


16 


0.7045 


96 


0.22911(72) 


0.1054(34) 


e-5 


-0.000019(11) 


lo.5 


0.73 


96 


0.2023(11) 


0.42(2) 


e-5 


0.000020(23) 


17 


0.7555 


64 


0.179(1) 


0.374(12) 


e-3 


-0.000022(82) 


17 


0.7555 


96 


0.17942(59) 


0.1393(37) 


e-4 


-0.000018(25) 


17.5 


0.7809 


96 


0.157(1) 


0.455(21) 


e-4 


0.0000057(46) 


18 


0.8063 


64 


0.13718(87) 


0.1783(49) 


e-2 


-0.00005(17) 


18 


0.8063 


96 


0.13898(63) 


0.1218(35) 


e-3 


-0.000006(48) 


18.5 


0.8317 


96 


0.11931(87) 


0.356(14) 


e-3 


-0.000027(92) 


19 


0.857 


64 


0.1005(10) 


0.761(26) 


e-2 


0.00073(52) 


19 


0.857 


96 


0.10025(46) 


0.1023(23) 


e-2 


-0.00005(12) 


19.5 


0.8824 


96 


0.08206(66) 


0.2895(94) 


e-2 


0.00010(33) 


20 


0.9077 


84 


0.0635(11) 


0.1476(86) 


e-1 


0.0017(18) 


20 


0.9077 


96 


0.0642(13) 


0.846(60) 


e-2 


0.0033(15) 


20 


0.9077 


160 


0.06332(7) 


0.699(37) 


e-3 


0.00015(24) 



Table 1: The values of the parameters in equation ([8]) fitted to the data points 
shown in the figure [TJ 
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T/ctq value shown, the data points exhibit an Ns dependence, to be associated 
with effects of the nearby critical temperature (see fllSp ). 

At this point, it must be emphasized that this critical temperature is aston- 
ishingly close to the location of the Nambu-Goto singularity, which is from Eq. 



As a last remark , we point out that while the string model is often supposed 
to be relevant to SU{N), iV — >• oo [8], the transition is first order for > 5 [21 [7j , 
and thus cannot be associated with the Hagedorn temperature of the string model. 

5 Conclusions 

In this article we have described a numerical calculation on the lattice of the 
correlations between Polyakov loops in three dimensional SU(3) gauge theory at 
finite temperature. We have chosen A^,- = 6, which we believe is large enough to be 
in the scaling region in the three dimensional theory. We vary the temperature 
by varying the lattice coupling constant. By doing this we get a dense set of 
measurements between T = Tc/2 and Tc. We also calculated the susceptibility 
of the Polyakov loops in the neighboorhood of the deconfinement transition, and 
the moments needed for the Binder method applied to the variable From 
this we obtain a determination of the critical lattice coupling constant. Using 
results from Teper et al. we construct an interpolating function and determine 
the dimensionless physical quantity T^j ^fa^ = 0.976(15). It is compatible with 
the Hagedorn temperature of the Nambu-Goto string model. This is somewhat 
astonishing, because the string model has the mean field exponent near the phase 
transition u = 1/2, whereas from universality arguments we expect u = 5/6. 
To investigate this further, we use the correlation lengths from the correlation 
between Polyakov loops at temperatures below the phase transition, where there 
is confinement and the string model may be applicable. The correlation lengths 
are inversely proportional to the string tension within the string interpretation. 
We find a very good agreement with the string tension from the Nambu-Goto 
string model up to T ^ 0.7Tc. However, up to this temperature the universal 
terms for a bosonic string, which coincide with the expansion of the Nambu- 
Goto model also describe the data. Above this temperature we see a systematic 
deviation from the string behaviour. This is due to the fact that here we enter 
the critical region. 

In a subsequent paper [12] we will study in detail the critical region of the decon- 
finement transition in the three dimensional SU(3) gauge theory. 



dZOD 




(23) 
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